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Abstract 

Recently, a number of authors have investigated the conditions under 
which a stochastic perturbation acting on an infinite dimensional dynami- 
cal system, e.g. a partial diflterential equation, makes the system ergodic and 
mixing. In particular, one is interested in finding minimal and physically nat- 
ural conditions on the nature of the stochastic perturbation. I shall review 
recent results on this question; in particular, I shall discuss the Navier-Stokes 
equation on a two dimensional torus with a random force which is white noise 
in time, and excites only a finite number of modes. The number of excited 
modes depends on the viscosity v, and grows like when v goes to zero. 
This Markov process has a unique invariant measure and is exponentially 
mixing in time. 
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1. Introduction 

The goal of this paper is to consider stochastic partial differential equations and 
to study conditions on the random perturbation that imply exponential convergence 
to a stationary state. In fact, one wants 'minimal' conditions, in the following sense: 
by expanding the solution in a basis of eigenfunctions of a linear operator associated 
with the PDE, one can write the latter as an infinite dimensional system of coupled 
differential equations. The question, then, is: to how many such equations do we 
need to add noise in order to make the system ergodic and mixing? 
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The physical motivation for this question comes from the fact that isotropic 
turbulence is often mathematically modelled by Navier Stokes equation subjected 
to an external stochastic driving force which is stationary in space and time. If the 
solution is expanded into Fourier modes, the driving force, which, in the language 
of physicists, acts on "large scale", should not perturb, or perturb very weakly, the 
high modes which represent the small scale properties of the system. So, one would 
like to show that the system becomes ergodic and mixing by adding noise to as 
few modes as possible. Obviously, this requires some detailed understanding of the 
nonlinear dynamics of the deterministic PDE. 

This problem is interesting from another point of view. As we shall see below, 
one can show that all but a finite number of modes converge to equilibrium provided 
the remaining ones do. So, we can reduce ourselves to a finite dimensional problem, 
which would be standard, except for the fact that the discarded modes produce a 
memory effect on the remaining ones, so that the problem is no longer Markovian. 
At this point, one introduces techniques coming from the study of the statistical 
mechanics of one-dimensional systems (where the unique dimension corresponds 
physically to space rather than time) with "long range, exponentially decaying, 
interactions" which have already been very useful in the study of SRB measures in 
dynamical systems (see p8||). 

At present, the best results require that the number of modes to which noise 
must be added depends on the parameters of the system, although a stronger result 
is likely to hold (see Remark 4 after Theorem 1.1). 

The type of question discussed here (for the Navier-Stokes equation but also 
for other equations) has been at the center of attention of several groups of people 
(see Remark 3 after Theorem 1.1 below). In this paper, I shall try to explain, in a 
simplified form, the approach followed by A. Kupiainen, R. Lefevere and myself in 
1^ (see also ||l|, ^ for previous results). 

To be concrete, consider the stochastic Navier-Stokes equation for the velocity 
field u(t,x) G defined on the torus T = (R/27rZ)2: 

du+ {{u-V)u- iyV^u + Vp)dt = df (1.1) 



where f(t,x) is a Wiener process with covariance 



EU{t,x)fp{t',y) = mm{t,t'}Ca,0ix - y) (1.2) 



and Cap is a smooth function satisfying daCap = 0. Equation (|l . l| ) is supple- 
mented with the incompressibility condition V • u = = V • /, and we will also 
assume that the averages over the torus vanish: w(0, x) = Q = J!p /(t, x), which 
imply that Jrp uit, x) = Q for all times t. 



It is convenient to change (1.1) to dimcnsionlcss variables so that v becomes 



equal to one. This is achieved by setting u{t, x) = vu'{vt, x). Then u' satisfies (1.1) 



(1.2) with u replaced by 1, and C by 
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From now on, we work with such variables and drop the primes. The dimcnsionless 
control parameter in the problem is the (rescaled) energy injection rate i-trC"(0) , 
customarily written as (Re)'^ where Re is the Reynolds number: 

Re = e^iy'^, 

and e = 5-trC(0) is the energy injection rate in the original units (for explanations 
of the terminology see pO|). 

In two dimensions, the incompressibility condition can be conveniently solved 
by expressing the velocity field in terms of the vorticity uj — diU2 — d2Ui. First 
(1.1) implies the transport equation 

duj + {{u ■ V)lu - V'^uj)dt = db, (1.3) 

where b — 81/2 — 82/1 has the covariance 

Eb{t, x)b{t',y) = mm{t, t'}(2Tr)-^-f{x ~ y) 

with 7 = -27r!/-3AtrC. 

Next, going to the Fourier transform, ujk(t) — Jrj, e^'''^uj{t,x)dx, with k G 
Z^; we may express u as — i-^— ^^'j-'fc, and write the vorticity equation as 

duj{t) = F{uj{t))dt + db{t), (1.4) 

where the drift is given by 

F(u;)fe = -fcV- + ^ ^-^^^^wu^iuji (1.5) 

iez2\{o,fc} 

and {&fe} are Brownian motions with b^ — and 

Ebk{t)bi{t') ^ TmTi{t,t'}5k,-nk. 
The dimcnsionless control parameter for the vorticity equation is 

^ = II 7fc = 2^7(0) (1.6) 

which is proportional to the ld injection rate, and also to the third power of the 
Reynolds number. One is interested in the turbulent region, where R is large; 
therefore, we will always assume below, when it is convenient, that R is sufficiently 
large. 

For turbulence one is interested in the properties of stationary state of the 
stochastic equation (L4) in the case of smooth forcing (see [0 for some discussion of 
this issue) and, ideally, one would like to consider the case where one excites only 
a finite number of modes, 

7fe ^ , < TV, 
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with N of order of one (for that, see Remark 4 below). In this paper we assume 
that N scales as 

N = kR, (1.7) 

with K a constant, taken large enough. We set all the other jk = 0, although this 
condition can easily be relaxed. Let us denote the minimum of the covariance by 

p-min{|7fe| I \k\''<N}. 

Before stating our result, we need some definitions. Let P be the orthogonal 
projection in iJ = L'^(T) to the subspace Hg of functions having zero Fourier 
components for |/cp > N. We will write 

UJ = s + I 

with s = Poj, I = (1 — P)uj (respectively, the small k and large k parts of oj). Denote 
also by Hi the complementary subspace (containing the nonzero components of ^). 
H is our probability space, equipped with the Borel cr-algebra. 

The stochastic equation ( |l.4| ) gives rise to a Markov process U!{t) and we denote 
by P*{E\uj) the transition probability of this process. 

The main result of 01 is the 



Theorem 1.1. The stochastic Navier-Stokes equation (1.4) defines a Markov 
process with state space {H,B) and for all R < oo, p > it has a unique invariant 
measure fi there. Moreover, Vo; G H , for all Borel sets E G and for all bounded 
Holder continuous functions F on Hi, we have, 

P\duj'\oj)lE{s')F{l') - I p{duj')lE{s')F{l'))\ < Cll^^lUe-™* (1.8) 



where C = C(||a;|| , i?, p) < oo, m = m{R,p,a) > 0, and Hi^Ha is the Holder norm 
of exponent a. 

Remark 1. Infl, we stated, for convenience. Theorem 1.1 by saying that the 



constant C in (L8) was a function of lu which was almost surely finite. Since this 



was stressed e.g. in [20 , it is worth remarking that C is simply a function of 
(depending also on the parameters R and p), which is finite Vw G H. To check this, 
we refer the reader to equations (86) and (97) in The main reason why this 
bound holds, however, lies in the fact that the only dependence of our estimates on 
Lo appears in Lemma 4.1 below and occurs through \\uj\\. 

Remark 2. In |l| it was proven that, with probability 1, the functions on the sup- 
port of the measure constructed here are real analytic. In particular all correlation 
functions of the form 

/i(dw) J|V"'u(a;,) 

i 

exist. For further results on analyticity, see E 
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Remark 3. While the existence of the invariant measure follows with soft methods 
||29| , its uniqueness and the ergodic and mixing properties of the process has been 
harder to establish. With a nonsmooth forcing (meaning that the strength of the 
noise, 7fc, decays only polynomially with \k\) this was established in and for 
large viscosity in p3| . However, those results did not cover the most physically 
interesting situations. The first result for a smooth forcing was by Kuksin and 
Shirikyan who considered a periodically kicked system with bounded kicks (for 
results on exponential convergence in that model, see |l^, |l^, ^ ) . In particular 
they could deal with the case where only a finite number of modes are excited by the 
noise (the number of modes depends both on the viscosity and the size of the kicks). 
In 1^ , we proved uniqueness and exponential mixing for such a kicked system where 
the kicks have a Gaussian distribution, but we required that there be a nonzero 
noise for each mode. An essential ingredient in analysis of jlSj, which was used 
in and by other authors, is the Lyapunov-Schmidt type reduction that allows 
to transform the original Markov process with infinite dimensional state space to 
a non-Markovian process with finite dimensional state space. While the analysis 
of [l^ was limited to bounded noise acting at discrete times, it was extended in 
jlfil, |l7|, [20I , I21I to cover unbounded noise and continuous time, as well as to obtain 
results on the strong law of large number and the central limit theorem. The first 
results on ergodicity of the system with unbounded noise and finitely many excited 
modes were obtained in |^ (see also for applications to other equations) and, 
for exponential convergence, in J^, which was also proved in p^ . For results on 
related problems, see |ll|, |l^, |2J. 

Remark 4. What one would like to obtain is a result similar to Theorem 1.1, 
but with TV finite, independently of R. An interesting result in that direction was 
obtained by Weinan E and Mattingly j|] who showed that, if one adds noise to only 



2 (suitably chosen) modes, ergodicity holds, provided one truncates the system ( |1.4 , 
|l.5| ), by keeping only a finite, but arbitrarily large, number of modes. This of course 
suggests that the 2 stochatically perturbed modes produce an "effective noise" on 
any finite number of modes, in particular on all those with < kR; then, one could 
hope to combine this with the results in |Q, |^ to obtain ergodicity and mixing for 
the full system. This, however, has not been done. 

Remark 5. The parameters in our problem are R and p. All constants that do not 
depend on them will be generically denoted by C or c. These constants can vary 
from place to place. 



Let me now explain the connection with ideas coming from statistical mechan- 



ics. 



First, observe that, if one neglects the nonlinear term in (1.4-1~5|), one expects 



Ijo^ll to be of order R^ , for typical realizations of the noise (i?^ is the typical size of 



the noise, and the —k-Lo^ term will dominate in eq. (1.4) for larger values of 



It turns out that similar probabilistic estimates hold for the full equation (1.4) as 



shown in Section 4. Now, if is of size R"^ , the —k'^ivk term will dominate the 



nonlinear term (which is roughly of size UtiJll ) in eq. ([l.4|), for |fc| > k_R^, and one 
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can expect that those modes (corresponding to / above) wiU behave somewhat Uke 
the solution of the heat equation and, in particular, that they will converge to a 
stationary state. 

Thus, the first step is to express the ^-modes in terms of the s-modes at 
previous times. This is done in Section 2 and produces a process for the s-modes 
that is no longer Markovian but has an infinite memory. In statistical mechanics, 
this would correspond to a system of unbounded spins (the s-modes) with infinite 
range interactions, with the added complications that, here, the measure is not given 
in a Gibbsian form, but only through a Girsanov formula, i.e. ( |2.9| ) below, and that 
time is continuous. Hence, we have to solve several problems: the possibility that 
uj be atypically large, the long range "interactions", and finally, showing that a 
version of the s-process with a suitable cutoff is ergodic and mixing. 

In Section 3, 1 introduce a "toy model" , namely a process with infinite memory, 
but with bounded variables, so that the problems caused by the unprobably large 
values of \\uj\\ does not occur. In that model, I explain how the statistical mechanical 
techniques, developed to study systems on one dimensional lattices, can be adapted 
to our setting. 

The large uj problem is treated in Section 4, using probabilistic estimates 
developed in which, in statistical mechanics, would be called stability estimates. 
In Section 5, I sketch how the remaining problems are handled: showing that the 
techniques explained in Section 3 can be applied here. However, this is where 
several technical complications enter, for the treatment of which I refer to The 
problem is that, even though for typical noise, hence for typical cj's, the Z-modes 
depend exponentially weakly on their past (see Section 2), thus producing, typically, 
"interactions" that decay exponentially fast, they may depend sensitively on their 
past when the noise is large. In the language of statistical mechanics, atypically 
large noise produces long range correlations, and that is the source of many technical 
difficulties. My goal here is to present the main conceptual tools used in g, putting 
aside those difficulties. 



2. Finite dimensional reduction 

Using an idea of one can reduce the problem of the study of a Markov 
process with infinite dimensional state space to that of a non-Markovian process 
with finite dimensional state space. 

For this purpose, write the equation ( |l.4| ) for the small and large components 
of w separately: 

ds{t) = PF{s{t) + l{t))dt + dh{t) (2.1) 
j-^m = {l-P)F{s{t) + l{t)). (2.2) 

The idea of is to solve the / equation for a given function s, thereby defining 
l{t) as a function of the entire history of s{t'), t' < t. Then, the s equation will 
have a drift with memory. Let us fix some notation. For a time interval /, we 
denote the restriction of ui (or s, / respectively) to / by lu{I), and use the boldface 
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notation s(/), to contrast it with s{t), the value of s at a single time. || • || will 
denote the norm. In [Q it was proven that, for any t < oo, there exists a set Br 
of Brownian paths 6 e C{[Q,t], Hs) of full measure such that, for b G Br, (1.4) has 
a unique solution with < oo, ||Vu;(t)|| < oo for all t (actually, u!{t) is real 

analytic). In particular, the projections s and I of this solution are in C([0, r], -ffs(i)) 
respectively. 

On the other hand, let us denote, given any s G C{[0,t], Hg), the solution 
— whose existence will be discussed below — of (2.2), with initial condition 1(0) 
by l{t,s{[0,t]),l{0)). More generally, given initial data l{t') at time t' < t and 
s([i',T]), the solution of (2.2) is denoted, for cr < r, by l{a,s{[t' ,a]),l{t')) and the 
corresponding oj by uj{(T,s{[t' ,a]),l{t')). The existence and key properties of those 
functions are given by: 

Proposition 2.1. Let 1{Q) e Hi and s e C{[0,t],Hs) . Then l{-,s{[0,t]),l{0)) e 
C{[0,t], Hi)nL'^{[0,T], Hi), where Hi ^ HiHH'^ , and H^ is the first Sobolev space. 
Moreover, given two initial conditions li,l2 and t <t 



\l{t, s([0, t]), h) - l{t, s([0, t]), l2)\\ < exp 



-KRt 



11^1-^211 (2.3) 



where a = {2n) ^X^l^l ^ 
satisfies 



and uji(t) = s{t) + ^i(i, s([0, t]), Zi). The solution also 
lit, s([0, t]), m) = lit, s([r, t]), liT, s([0, r]), Z(0))). (2.4) 



Remark. What this Proposition shows is that the dependence of the function 
I upon its initial condition li, i = 1,2, decays exponentially in time (i.e. like 
the solution of the heat equation), provided uj is not too large, in the sense that 
J* IjVwilP < cRt, for a suitable constant c. As we will see in Section 4, this event 
is highly probable. 

Now, if s = Plu with uj as above being the solution of ( |1.4| ) with noise b € Br 
then the /(s) constructed in the Proposition equals (1 — P)llj and the stochastic 
process s(i) satisfies the reduced equation 

dsit) ^ fit)dt + dbit) (2.5) 

with 

fit) = PFicuit)). (2.6) 
where w(t) is the function on Ci[Q,t], Hg) x Hi given by 

u;it) = sit) + lit,si[0,t]),liO)). (2.7) 



(2.5) has almost surely bounded paths and we have a Girsanov representation for 



the transition probability of the w-process in terms of the s-variables 

P*(F|^(0))- / ML(o)(ds)i^(^(t)) (2.; 
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with 

^^Uo^ds) = e/o(/(^)--H.s(.)-i/M..)),*^^^(rf,) (2.9) 

where i^I^^q^ is the Wiener measure with covariance 7 on paths s = s([0, i]) with 
starting point s(0) and (•, •) the £^ scalar product. Define the operator 7"^ in terms 
of its action on the Fourier coefficients: 

if,i-'f)= E (2-10) 

|fc|2<JV 

The Girsanov representation (|]^) is convenient since the problem of a stochas- 
tic PDE has been reduced to that of a stochastic process with finite dimensional 
state space. The drawback is that this process has infinite memory. In the next 
section, I will show how to deal with this problem in a simplified situation. 



3. A Toy Model 

In order to explain the main ideas in the proof, I will consider first a 'toy 
model' and then explain the steps needed to control the full model. 

Let us consider variables Xt G [0, 1], t € Z about which a set of (consistent) 
conditional probability densities p{xt\:x.[^oc,t-i]) is given, i.e. one is given the prob- 
ability densities of the variables Xt, at time t, given a 'past history' ^[-oo,t-i]^ where 
we write, for / C Z, x/ = {xt)tei & [0, 1]"'^. 

Before stating precise assumptions on p, here is what one wants to prove: 
3C < 00, m > and a probability p on [0, 1] such that VE C [0, 1], E measurable, 

\p{xT e E\^[_^,o])~P{E)\ < Ce-"^ (3.1) 

for all T > and all x[„oo,o]i where 

^1 T-l T 



•^0 t=i t=i 



..0]) (3.2) 



and X[i V X[_oo 0] denotes the obvious configuration on [—0 0, t — 1]. 

Now let us state the assumptions on p that will imply (3J); obviously, we 
assume that: 

p(a;t|x[_oo,t-i]) > (3.3) 

and 



L 



1 

dxtp{Xt\X[-ac>,t-l]) = 1 (3.4) 



for all X[-oo,t-i]- Moreover, we assume that p{-\-) is invariant under translations of 
the lattice Z, in a natural way. The non-trivial assumptions are: 
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a) Let, for s < t — 2, 

Ss,t{:>^[s,t]) = P{xt\^[s,t-i] VO) - p{xt\x[s+i^t-i] VO), (3.5) 

where x/ V denotes the configuration equal to xt for t I and equal to zero 
elsewhere. We assume that 3C < 0, m > such that Vs, i e Z as above, 

llf^s.tlloo < Cexp(-m|t-s|). (3.6) 

b) Define, for iV > 1, the Markov chain on fl = [0, 1]^ by the transition 
probability 

N 

^'(x[i^Ar]|x[_Ar+i,o]) = J]^ J3(a;t |x[t_jv^t„i] V 0). (3.7) 

We assume that this Markov chain satisfies : 3S > 0,Vi? C ri,Vx, x' G fl, 

P(B|x) +P(S"|x') > 5 (3.8) 

where S is independent of N (see however the Remark following the proof of Propo- 
sition 3.1 for a generalization). 



Proposition 3.1. Under assumptions a) and h) above, (3.1) holds. 



Remark 1. The techniques used here can also prove the analogue of ( ^.l| ) with xt 
replaced by 2^[t,t-l]i for ^^ly finite L, and this, in turn, allows one to associate to 
the system of conditional probabilities a unique probability distribution on [0, 1]^ 
(which is called, in statistical mechanics, the Gibbs state associated to the system 
of conditional probabilities), but I will not go into that, because I want to give here 
only an elementary idea of the techniques used in Of course, this type of results 
is not new (see e.g. ||2^, Lecture 12, for a similar result, applied to dynamical 
systems, with a somewhat different proof). 

To prove the Proposition, we first use a resuh of Doob (g, p. 197-198): 

Lemma 3.1 For the Markov chain defined in h) above, there exists a probability 
distribution P on Q such that Vx G Q, V-B C fi, Vn > 1, 

|P"(B|x)-P(S)| < (l-<5)". (3.9) 

Proof. Let P{B,n) = supxP"(P|x) and P{B,n) = infxP"(P|x). It is easy to 
see that P{B,n) is decreasing in n, while P_{B,n) is increasing in n. Thus, it is 



sufficient to prove the bound (3.9) for the difference |P(P,n) — P.{B,n)\ and, for 
that, we shall prove: 

< P(P, n + 1) - P(P, n + 1) < (1 - 5)(P{B, n) - P{B, n)). (3.10) 



Since P{B,n) - P{B,n) < 1, follows. 
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Define a signed measure on subsets of ft: 

*x,x'(^^) - PiE\x.) - P{E\x') (3.11) 

and let S+ (resp. S^) denote the set where 5'x,x'(£') > fov E C S+ (resp. < 0). 
We have: 

P{B,n+ 1) - P{B, 71+ 1) ^ sup [ [P(dx"|x)-P(dx"|x')]^"(S|x") 

x,x' J 

- sup /'M'x,x'(dx")P"(S|x") 

x,x' J 

< sup(4'x,x'(5'+)P(P,n) + *x,x'(5-)P(B,n)). (3.12) 

x.x' 

By definition, ^'x,x'(S'") = -^'x,x'('S'^): so tliat 

*x.x'(5+)P(P, n) + *x.x'(^")P(P, n) = *x,x'(5+)(P(P, n) - P{B, n)). 
Also, for any set P C il, ( |3.8[ ) implies 

*x,x'(£;) = 1 - {P{E'\^) + P(P|x')) < 1 - 
Applying this to P = 5*+ in (|1|) implies ( |3lo| ). 



Remark 2. We shall use this Lemma under the following form: 

[ dx|P"(x|x')-P(x)| <2(l-(5)" (3.13) 
Jn 



for all x' S ri; this follows by applying (3^) separately to the sets where the 
integrand is positive and negative. 
Now, let us turn to the 

Proof of Proposition 3.1. 

We write each factor in (13.21) as 



p(xt|x[_oo,i-i]) =p(a;t|x[t_Ar,t_i] V0)+ ^ (5s,t(x[s,t]), (3.14) 

\s-t\>N 



where N is an integer to be chosen later. Insert this in the product in (3.2), and 
expand: we get 

T-l 

P(xt|X[_oo,0]) = 

Yl Yl / \{'i^t\{^s,M[sA)^P^^t\'^[t-N.t-i]yQ), (3.15) 
/c[i,T] s •' t=i tei t^i 



where the sum over subsets / corresponds to the choice in (|3.14| ) between the first 
term and the sum, while the sum over s = (st)tg/ corresponds to the possible choices 
of a term in that sum. 
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Now, let / = (J [st , t] and let 
tei 

[l,r]\7-|Jj,|Jj„, (3.16) 

i a 

where each Ji is a union of intervals of length N, containing at least two such 
intervals, and each la is an interval of length less than 2N between two connected 
intervals in / or an interval of length less than N between an interval in / and 



an interval Ji. The reason for these definitions is that, in the RHS of ( 3.15 ), the 
only functions depending on Xs, with s in the complement of /, arc the factors 
p{xt\x[t-N .t-i]) J so that, by integrating over these variables, one can obtain the 
transition probabilities of the Markov chain defined in condition b) above. For 
that, we need intervals of length at least 2N, which are the J^'s, while the intervals 
/q,'s simply cover the leftover sites. 

Since the model here is translation invariant, let us fix one interval Ji — J, 

n 

and write it as a union of disjoint intervals of length N: J — \^Ki with Ki = 

1=0 

[t + l + lN,t+{l + l)N]. 

We have, by definition (3.7) of the transition probability P: 

/n— 1 n 
n ^^A-' n n {p{^M[t-N.t-i] v o)) = p^^kj^ko)- (3.17) 
1=1 l = lt£Ki 

Now write this as 



P"(xkJx;^J - P(xA'J + P(xkJ, 



(3.18) 



where P is defined by (B.£). Apply this to each interval Ji in (3.16), with n re- 
placed by rii — — 1. Insert that identity in ( 3.15 ) for each Ji and expand the 
corresponding product over i oi Ai + Bi, where A = P"(x/f^|x/f„) — P{xk„) and 
P = P(xkJ. 

For E as in (|]^), integrate over E each term in the resulting expansion, and 

write 



p(xTeP|x[_oc,o]) = Q + ^, (3.19) 

where Q collects all the terms in the resulting sum where at least one factor P{xk^. ) 
appears and R all the rest. Now, the presence of one such factor P 'decouples' xt 
from the initial conditions x^-oofl], in the sense that, if we consider the difference 

pixT e P|x[_^,o]) --pixT e £:|x{_^ o]), (3.20) 

for two different past histories, then the Q sums are equal and only the R sums 
contribute to the difference. Indeed, fix a and consider all the terms in our 
expansions where the factor P(x/<-^ ) appears; let tg be the last time before the 
interval _?C„ . . By construction, in all the terms under consideration, all the functions 
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that depend on xt, for t > to do not depend on the variables Xt, for t < tQ. So, if we 
resum, in the expansion, all the terms depending on the variables Xt , for t <tQ, we 
obtain, for the two terms in ( |3.20 ), p(a;to |x[_oo,o]): ^-iid |x'_^ ^j) (we simply 
use (3.15) read from right to left, with T replaced by to). But performing in (3.2) 
the integral over xtg gives 1 in both cases, which shows that the difference between 
the respective sums cancel. 

So, if we show that, 3C < oo,m > such that 



\R\ < Ce-™^, 



(3.21) 



Vx[_Qo,o]j we obtain that the absolute value of (3.20) is exponentially small and, 
from that, (3.1) easily follows. 



(3.22) 



Using the bound ( |3.6[ ) on 6s,f and (3.13) on 

/ rr dxt\p'\xK„\^Ko) - p{^K„ 



and the fact that, by (3.4) and xt £ [0, 1], all the integrals are bounded by 1, , we 
get: 



1^1 ^ E E n(^^^™'*~''' ) 11(2(1 - 6r) 



(3.23) 



s tei 



where the second product runs over the intervals Ji in (3.16), and where rii — — 1. 
Note that the length of each la in ( |3.16| ) is less than 2N and, since such intervals 
are always adjacent to a connected component of / (unless / = 0, in which case 
this number is at most 2), the number of inter vals I g is less than 2|/| + 2; the same 
bound holds for the number of intervals Jt in ( 3.16 ) (in fact, a better bomid holds 
here, but we won't use it). So, we have: 



N 



2) > 



(T- 



N 



c\I\ 



(3.24) 



for some number c, where, in the second inequality, we use < 2N and ( [sl^ ). 
Using this, we can, by changing the constant C, bound ( [3.23| ) by: 



^ E E n e-™l*-^'lcl^l (1 - 5)(^-l^ 



\)/N 



(3.25) 



Since by definition of /, ^~^|^ ^ St\ > |-^|, we can, by considering separately the 



terms where |/| < -j, and those where |/| > ^, bound the sum in (3.25) by 



where 



m = mm 



m — ln(l — S) 
T' 2N 



(3.27) 
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Now, choose N so that 



\t-s\>N 



with 



(I + C77) < e 



m/2 



(3.28) 



(3.29) 



which is po ssible since, from ( 3.28 ) we see that, for large N, rj = exp{—0{N)) while, 
from (^!27|) , m = 0{N-^). 

We use ( 3.28| ) to control the sum over each st in (|3.26| ), and we get 



( p6| ) < Ce"*^ J2 (<^^)'^' < C'(l + Crj) 

/C[1,T] 



T^-rnT 



(3.30) 



and, using (3.2£), we get ( 3.2l| ) with 



This completes the proof of Proposition 3.1. 



Remark 3. By considering ( ^^2^ , |3.28| , ^.29| ), we sec that one can extend the proof 
to a situation where ^ in ( 3.S ) depends on N, as long as (5 > exp(— cA^) for a constant 
c small enough. 

Now, let us turn to the real model, and make a list of the difficulties not 
present in the toy model. The first one is that time is continuous rather than 
discrete, but that is a minor problem. We can easily introduce a discretization of 
time. A more serious problem is that one deals with what are called "unbounded 
spins" in statistical mechani cs o r what is also known as a "large field problem", 
namely the variables s(r) in ( ^.9[ ), which play a role similar to the variables Xt here, 
take value in rather than [0, 1] (actually, if we consider the variable s over a 
unit time interval, they take values in a space of functions from that interval into 
R^). And, what really causes a problem, is the fact that the bounds (3^), ( |3.8[ ) 
do not hold when the variables s take large values. However, as we sha ll see in th e 
next section, this is unprobable. Thus, before doing an expansion as in ( ^.15| , 3.18| ), 
we must first distinguish between time intervals where the s variables are large and 
those where they are small. Then, putting aside lots of technicalities, we perform 
the expansion (3.15) in the latter intervals and use estimates like (4.2) below to 
control the sum over the intervals where co is large. 

Finally, there is an additional difficulty coming from the fact that the definition 
of the probabilities here involve a Girsanov representation. In statistical mechanics, 
one usually deals with situations where the probabilities (^.3D can be written as: 



p{xt\x[^aa,t~i]) = exp(E (j)i{xi)), 



(3.31) 



where the {0/}'s represent "many body interactions' (suitably normalized so that 
(3.4) holds) and the sum runs over intervals / C Z whose last point is t. Then, a 
bound of the form 



\(j,i\\oo < Ccxp(-m|/|), 



(3.32) 
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with C < oo, m > 0, is enough to obtain ( |3.6D and ( |3.8| ). But here the probabihties 



are not of that form, because of the stochastic integral J^^ /(t)7 ^c?s(r) in (2.9) 



4. A priori estimates on the transition probabili- 
ties 



The memory in the process ( 2.5) is coming from the dependence of the solution 
of ( |2.2[ ) on its initial conditions. By Proposition 2.1, the dependence is weak if 
/p llVwp is less than cRt for a suitable c. It is convenient to define, for each unit 
interval [n — 1, n] = n, a quantity measuring the size of w on that interval by: 

sup Mt)\f + f \\Wiu{t)\fdt. (4.1) 

ten Jn 

The following Proposition bounds the probability of the unlikely event that 
we are interested in: 

Proposition 4.1. There exist constants c > 0, c' < oo, (3q < oo, such that for all 
t,t',l<t < t' and all 13 > (3o, 

t'-i 

P{J2 Dn(uj) > PR\t' - i||c^(0)) < exp(lc'e-*|lu;(0)|p) exp{-c(3\t' ~ t\). (4.2) 



Remark 1. This means that the probability that w is large over an interval of time 
decays exponentially with the length of that interval, provided that ||w(0)|| is not 
too large. And, if ||w(0)|p is of order Dn(y>) will be, with large probability, of 
order R after a time of order logX. 

The main idea in the proof is a probabilistic analogue of the so-called enstrophy 
balance: in t he d eterministic case, using integration by parts and V • u = 0, on 
derives from (L3) with dh = 0, the identity: 

which implies that the enstrophy (||w||^) decreases in time. This basic property of 
equation (|l.3|) makes the proof of the following Lemma rather simple. 



Lemma 4.1. For all cli(O) G L"^ , and all t > 0, 



E 



< 3e' 



(4.3) 



and 



P{\\uj{t)\\^ > D\uj{0)) < 3e-3s, 



"'11^(0)11' 



(4.4) 



Ergodicity and Mixing for Stochastic Partial Differential Equations 



581 



Proof. Let x{t) = X{t)\\uj{t\P = KT)J2k\^k\'^ ^r < r < i. Then by Ito's 



formula (remember that, by ( |l.6| ), 7^ = R and thus 7^ < R, Vfc): 



A5]7fc + 2A2^7fekfene^ 



< E[{{X\-^ - 2 + 2\R)x + \R)e% 



(4.5) 



where E denotes the conditional expectation, given a;(0), and where we used the 
Navier-Stokes equation ( |1.3D , |fc| > 1 for ujk 7^ 0, and the fact that the nonlinear 
term does not contribute (using integration by parts and V • u = 0). Take now 
A(t) = A-e^^"*^ so that X<^, AA^i = 1, AA^^ - 2 + 2XR < -i and XR < \. So, 



i;[e^l < E\{^ - ia;)e^l < i - -Ele'V 
L J - LV^ 2 ' ^ - 2 4 ^ 



dr 



where the last inequality follows by using (1 — 2x)e^ < 2 — e^. Thus, Gronwall's 
inequality implies that: 

£;[e^(^)] < e-7e^(°) + 2 < 3e^(°), 

i.e., using the definition of A(t), 

^e-*||Lj(0)|i2 



E 



exp( — ||a;(r)|p) 



< 3 exp(- 



4i? 



)• 



This proves (4.3) by putting t = t; (4.4) follows from (4.3) by Chcbychev's inequal- 
ity. 



Since the _D„ in (4.2) is the supremum over unit time intervals of 



1 , 



||Vcj|pdT n~l<t<n, (4.6) 

In-l 

which does not involve only ||cj(t)|p, we need to control also the evolution of Dt{uj) 
over a unit time interval, taken, for now, to be [0,1]. From the Navier-Stokes 
equation (O) and Ito's formula, we obtain 



Dticj) = Doiuj) + Rt + / {uj,db) 



(since the nonlinear term does not contribute, as in (4.5)). 
Our basic estimate is: 

Lemma 4.2. There exist C < 00, c > such that, > 3Do{lj) 



P( sup Dt{cj) > A\ijiO)) < Cexp(-^). 
*e[o,i] 



(4.7) 



(4.8) 
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Remark 2. While the previous Lemma showed that tends to decrease as 

long as it is larger than 0{R), this Lemma shows that, in a unit interval, Dt{uj) does 
not increase too much relative to Dq{uj) = ^11^(0)11^. Thus, by combining these two 
Lemmas, we see that Dn{i^) — sup Dt{uj) is, with large probability, less than 

tG[n-l,n] 

||aj(0)|p, when the latter is larger than 0{R), at least for n > tiq not too small. 
Thus, it is unlikely that Dn{uj) remains much larger than R over some interval of 
(integer) times, and this is the basis of the proof of Proposition 4.1. 



Without entering into details, here are the main ideas in the proof of (4.8) 



From (4.7), we see that it is enough to get an upper bound on 



P sup I / {uj,db)\>{A- Do- R) uj{0) ] . (4.9) 

\tG[04] Jo / 

We use (see |^ for more details) Doob's inequality (see e.g.j27|, p. 24), to reduce the 
control over the supremum over t to estimates on | /q^(w, db)\. Letting E denote the 
conditional expectation, given w(0). and using Novikov's bound (see e.g. the proof 
of Lemma 5.2 below), we get 

/ .1 X 1/2 .1 s 1/2 

< i^j dTi;(e2^'("(")-^"M))j < i^J dr£;(e2^'^ll"MII')j (4.10) 

where the last two inequalities follow from Jensen's inequality, applied to 
^2e^ J^dr{u^{T),',uj(T))^ and from < R (see (O)). Now, choosing e so that 2e^R 



j^, i.e. £ = we can use (4^) to bound the RHS of ( 4.1C ). Combining this 



with Chebychev's inequality gives bounds on (|4.9| ) 



5. Decoupling estimates 



In this section, I shall give a very brief sketch of the ideas used to prove the 

of section 3 in the present setting, at least in 



analogue of assumptions (3.6) and (3 



the probable regions where co is small. The main point is to understand the analogue 
of the bound (3.32), which expresses the exponential decay of interactions. What 



plays the role of the right hand side of (3.31) is, see (2.9): 



9t 



(5.1) 



where, for simplicity, I consider a unit time interval [t — We want to show 

that this depends weakly on the past; so consider two functions gi, g2, defined 
in terms of t wo functions /i, /2, themselves defined through different li and I2 
(see (|2.6| , 2.7)). And, by analogy with what we did in section 3, we choose li = 
l{t, s([0, t]),l{0) = 0), I2 = l{t, s([l, t]),l{l) = 0), i.e. we set the large k modes equal 
to zero at different times (0 or 1). Using ( |2.4|) , we see that li — Z(i, s([l, t]), Zi(l)), 
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with — s([0, 1]), Z(0) — 0), so that we have, at time t = 1, two initial 

conditions, '2(1) = 0, with ||/i(l) — ^2(1)11 = ||'i(l)|| of order one, if ui is small 

in the interval [0, 1]. 

Now, if gt depends weakly on the past, it should mean that, for large t, gi and 
(72 are, in some sense, exponentially close. To measure the difference, write: 

ffi-52-(l--)gi = (l-i^)5i, (5.2) 
91 

which will be convenient, since we deal with unbounded variables for which sup 
norm estimates like in ( 3.32 ) are not available. Explicitly: 

jj = Jl,iSf{t),f-\dsit)-Mt)dt))~il*_^(Sf{t).f-^Sf{t))dt ^53^ 

where (5/ = /2 — /i- What we want to show is that 1 — iJ is, in a suitable sense, 
small. 

The next Lemma gives a bound on \\Sf\\ in terms of \\6l\\, and \\6l\\ is 

controlled by Proposition 2.1, provided that uj is small, in the sense discussed in 
section 4, in which case ||a;|j is also controlled, using sup^^^ 11^^(^)11 1^ (2I?n)^. 

Lemma 5.1. Let f{uj) — PF{uj) and uj = s + I, ui' = s + I' . Then, 

\\Sf\\ = \\f{^) - f{^')\\ < C(i?)(2||c^||||<5?|| + \\6ir) (5.4) 

with 51 = I — I' and C{R) a constant depending on the parameter R (see ( [j.6[ jj. 

Proof. We have 

^ Ifcl 

which, since < \J nR is bounded by 

VkR^^ \sk-pSlp + Sp6lk-p + Iplk-p - I'pl'k-pl- (5-5) 



Writing Iplk-p — lpl'^-p — IpSlk-p + lk-pSlp — Slpdlk-p and using Schwarz' inequality, 
we get 



(5.5) < VkR{2\\uj\\\\51\\\\ + \\Sl\\') 



which proves (5.4), since fk ^ only for k < kR, so that the sum in the norm 
1 1 (5/ 1 1 runs over C{R) terms. 

This Lemma would be enough to control (1 — H) if we had only in (5.3) the 
factor e^'^ Sf{t))dt^ which involves only ordinary integrals. 

To control the stochastic integral, it is convenient to undo the Girsanov trans- 
formation, i.e. to change variables from s back to b. Let E denote the expectation 
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with respect to the Brownian motion b with covariance 7 on the time interval [t—1, <]. 
We get, using (2.5): 

H = eA-iWW^T-'dbW)-^ fl,iSf(t),i-^5f(t))dt_ ^5 

Write now (1 — H)^ = 1 — 2H + to give a flavour of the estimates, let us 
see how one could show that the expectation with respect to E of —2H + H^is close 
to 1, i.e. that the expectation of (1 — i?)^ is close to zero. One can rather easily 
bound from below the expectation of H, using Jensen's inequality; to get an upper 
bound on the expectation of one uses: 

Lemma 5.2. Let C,{t) e C([0, l],-ffs) be progressively measurable. Then 

£:e/o(C.')'-'db)+A/„HC,7-k)<it < g2(l+A)||C||V-' (5 7) 

^;iere||C||=sup,||C(r)||2. 

Proof. This is just a Novikov bound: we bound the LHS, using Schwarz' inequality, 

by 

(^e^o(2C:7-'dfc)-2j(;(C,7-kdt))i(^g2(l + A)/oi(C,7-\)rft)i 

and note that the expression inside the first square root is the expectation of a 
martingale and equals one. 

We can then apply this Lemma to C = 2(5/, A = — i, replacing [0, 1] by [t— 1, <], 
and use the estimates coming from Lemma 5.1 and Proposition 2.1 to show that 
the RHS of (5/7) is exponentially close to 1, for i large. This gives a rough idea of 
why the "interactions" here are exponentially decaying, but it must be said that 
the full story is far more complicated and I refer to reader to for more details. 

Acknowledgments. 

I wish to thank my coauthors, Antti Kupiainen and Raphael Lefevere without 
whom this work would not have been possible. This work was supported in part by 
ESF/PRODYN. 



References 

[1] J. Bricmont, A. Kupiainen, R. Lefevere, Probabilistic estimates for the two di- 
mensional stochastic Navicr-Stokcs equations. J. Stat. Phys. 100 (3/4), (2000), 
743-756. 

[2] J. Bricmont, A. Kupiainen, R. Lefevere, Ergodicity of the 2D Navier-Stokes 
Equations with random forcing. Commun. Math. Phys. 224 (2001), 65-81. 

[3] J. Bricmont, A. Kupiainen, R. Lefevere, Exponential mixing of the 2D stochas- 
tic Navier-Stokes dynamics Commun. Math. Phys. (to appear). 

[4] J.L. Doob, Stochastic Processes^ John Wiley, 1953. 

[5] J. P. Eckmann, M. Hairer, Uniqueness of the invariant measure for a stochastic 
PDE driven by degenerate noise, Commun. Math. Phys. 219 (2001), 523-565. 

[6] W. E, D. Liu, Gibbsian dynamics and invariant measures for stochastic dissi- 
pative PDEs, J. Stat. Phys. (to appear). 



Ergodicity and Mixing for Stochastic Partial Differential Equations 



585 



[7] W. E, J.C. Mattingly, Ya.G. Sinai, Gibbsian dynamics and ergodicity for the 
stochastically forced Navier-Stokes equation, Comm. Math. Phys. 224 (2001), 
83-106. 

[8] W. E, J.C. Mattingly, Ergodicity for the Navier-Stokes equation with degen- 
erate random forcing, finite-dimensional approximation. Comm. on Pure and 
Appl. Math. LIV (2001), 1386-1402. 
[9] F. Flandoli, B. Maslowski, Ergodicity of the 2-D Navier-Stokes equation under 
random perturbations. Cow.mun. Math. Phys. 171 (1995), 119-141. 

[10] U. Frisch, Turbulence, Cambridge University Press, 1995. 

[11] M. Haircr, Exponential mixing properties of stochastic PDE's through asymp- 
totic coupling, Preprint (2001). 

[12] M. Haircr, Exponential mixing for a stochastic PDE driven by degenerate 
noise Nonlinearity 15 (2002), 271-279. 

[13] S.B. Kuksin, A. Shirikyan, Stochastic dissipative PDE's and Gibbs measures. 
Commun. Math. Phys. 213 (2000), 291-330. 

[14] S.B. Kuksin, A. Shirikyan, A coupling approach to randomly forced nonlinear 
PDE's. I, Comm. Math. Phys. 221 (2001), 351-366. 

[15] S.B. Kuksin, A. Piatnitski, A. Shirikyan, A coupling approach to randomly 
forced nonHnear PDE's. II, Comm. Math. Phys. (to appear). 

[16] S.B. Kuksin, A. Shirikyan, Ergodicity for the randomly forced 2D Navier- 
Stokes equations. Math. Phys. Anal. Geom. 4 (2001), 147-195. 

[17] S.B. Kuksin, A. Shirikyan, Coupling approach to white-forced nonlinear 
PDE's, Joum. de Math. Pures et Appl, 81 (2002), 567-602. 

[18] S.B. Kuksin, A. Shirikyan, On dissipative systems perturbed by bounded 
random kick- forces Ergodic Theory and Dynamical Systems (to appear). 

[19] S.B. Kuksin, On exponential convergence to a stationary measure of nonlinear 
PDEs, perturbed by random kick-forces, and the turbulence-limit. The M.I. 
Vishik Moscow PDE seminar, Amer. Math. Soc. Transl., (to appear). 

[20] S.B. Kuksin, Ergodic theorems for 2D statistical hydrodynamics, Rev. Math. 
Phys. (to appear). 

[21] S.B. Kuksin, A. Shirikyan, Some limiting properties of randomly forced 2D 
Navier-Stokes equations, preprint (2002). 

[22] N. Masmoudi, L.-S. Young, Ergodic theory of infinite dimensional systems with 
applications to dissipative parabolic PDE's. Preprint (2001). 

[23] J.C. Mattingly, Ergodicity of 2D Navier-Stokes equations with random forcing 
and large viscosity, Commun. Math. Phys. 206 (1999), 273-288. 

[24] J.C. Mattingly, Exponential convergence for the stochastically forced Navier- 
Stokes equations and other partially dissipative dynamics, preprint (2001). 

[25] J. C. Mattingly, The dissipative scale of the stochastic Navier-Stokes equation: 
rcgularization and analyticity, J. Stat. Phys (to appear). 

[26] A. Shirikyan, Analyticity of solutions for randomly forced two-dimensional 
Navier-Stokes equations, Russian Math. Surveys, (to appear). 

[27] B. Simon, Functional Integration and Quantum, Physics, Academic Press, 1979. 

[28] Ya.G. Sinai, Topics in Ergodic Theory., Princeton University Press, 1994. 

[29] M.I., Vishik, A.V. Fursikov, Mathematical Problems of Statistical Hydrodynam- 
ics, Kluwer, 1980. 



